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Abstract. The IQj symbol is a spin network that appears in the partition 
function for the Barrett-Crane model of Riemannian quantum gravity. Ele- 
mentary methods of calculating the lOj symbol require 0{j^) or more oper- 
ations and 0{j^) or more space, where j is the average spin. We present an 
algorithm that computes the lOj symbol using 0{j^) operations and 0{j^) 
space, and a variant that uses 0{j^) operations and a constant amount of 
space. An implementation has been made available on the web. 



1. Introduction 

The Barrett-Crane model of four-dimensional Riemannian quantum gravity 
has been of significant interest recently ||l|, ^, |l^, 0. The model is discrete and 
well-defined, and the partition function for the Perez-Rovelli version has been rigor- 
ously shown to converge Jill for a fixed triangulation of spacetime. The Riemannian 
model serves as a step along the way to understanding the less tractable but phys- 
ically more realistic Lorentzian version 0. However, despite its simplicity, we are 
currently lacking explicit numerical computations of the partition function and of 
expectation values of observables in the Riemannian model. These are necessary 
to test its large-scale behaviour and other physical properties. 

It has been shown |Q that the amplitudes in the Barrett-Crane model are al- 
ways non-negative, and therefore that the expectation values of observables can be 
approximated using the Metropolis algorithm. This greatly reduces the number of 
samples that must be taken, and thus the remaining obstacle is the time required to 
compute each sample. This paper presents a very efficient algorithm for doing these 
computations. The algorithm is used in Q and Q to understand the asymptotic 
behaviour of the lOj symbols and the dependence of the partition function on a 
cutoff. 

To explain further, we need to describe the Barrett-Crane model in more detail. 
It has been formulated by Baez 1^ as a discrete spin foam model, in which faces 
in the dual 2-skeleton of a fixed triangulation of spacetime are labeled by spins. 
The dual 2-skeleton consists of a dual vertex at the center of each 4-simplex of the 
triangulation, five dual edges incident to each dual vertex (one for each tetrahedron 
in the boundary of the 4-simplex), and ten dual faces incident to each dual vertex 
(one for each triangle in the boundary of the 4-simplex). 

Baez notes that the partition function for this model is the sum, over all labelings 
of the dual faces by spins, of an expression that contains the product of a lOj symbol 
for each dual vertex. A lOj symbol, described in detail in Section^, is a Spin(4) spin 
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network. Roughly speaking, a spin network is a graph whose vertices are labelled 
by tensors, and whose edges indicate how to contract these tensors. A spin network 
evaluates to a complex number in the way explained in Section ^ In short, the lOj 
symbol is a function taking ten input spins and producing a complex number. It is 
at the heart of the calculation of the partition function, and thus an algorithm for 
calculating the lOj symbols efficiently is quite important. 

In Section ^ we recall the definition of the lOj symbol using spin networks. 
Then in Section || we briefly describe the elementary algorithms for evaluating 
these spin networks, and give their running times and memory use. We conclude 
with Section ^, which presents our algorithms and their time and space needs. 

2. The lOj symbol 

In the dual 2-skeleton of a triangulation of a 4-manifold, each dual vertex belongs 
to five dual edges, and each pair of these dual edges borders a dual face. A lOj 
symbol is a Spin (4) spin network with five vertices (corresponding to the five dual 
edges) and ten edges, one connecting each pair of vertices (corresponding to the 
ten dual faces), with the edges labeled by spins. In the context of a Spin(4) spin 
network, a spin j labeling an edge denotes the representation j (g) j of Spin(4) = 
SU(2) X SU(2), where j is the spin-j representation of SU(2). Such representations 
are called "balanced." We use the convention that spins arc non-negative half- 
integers. 

Here is a picture of the lOj symbol, with the vertices numbered through 4, and 
the spins divided into two groups: are the spins on the edges joining vertex i 
to vertex i + 1 (modulo 5), and j2A are the spins on the edges joining vertex i to 
vertex i + 2 (modulo 5). 




The five vertices of the network are equal to Barrett-Crane intertwiners. These 
are the unique intertwiners (up to a factor) between four balanced representations 
of Spin(4) with the property that their expansion as a sum of tensor products of 
trivalent SU(2) networks only contains balanced representations on the internal 
edge, regardless of which pairs of external edges are joined Q. Barrett and Crane 
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give the formula for these intertwiners in 



Jl 






Here the sum is over ah admissible values of I, i.e. those that satisfy the Clebsch- 
Gordan condition for both SU(2) vertices. So I ranges from max(|ji — lis ~ jil) 
to min(ji + j2, js + ji) in integer steps. If the difference between these bounds is 
not an integer, the Spin(4) vertex will be zero. When I satisfies these conditions, 
there is a unique intertwiner up to normalization which can be used to label the 
trivalent SU(2) vertices. (These intertwiners are normalized so that the theta net- 
work in the numerator of equation (^) has value 1.) A; is the value of a loop in 
the spin-Z representation, which is just (— 1)^'(2Z + 1), the superdimension of the 
representation. 

The uniqueness result of Reisenberger tells us that if we replace the vertical 
edges in the above definition by horizontal edges, the result differs at most by a 
constant factor. In fact, Barrett and Crane |q| s tated that the two definitions give 
exactly the same Spin(4) vertex, and Yetter ]14| has proved this. 

Any closed spin network evaluates to a complex number, by contracting the 
tensors at the vertices according to the pairings specified by the edges. Thus the 
lOj symbol is a complex number. (In fact, one can show that it is always a real 
number.) 

To avoid confusion, we want to make it clear that we are working with the 
"classical" (non g-deformed) evaluation of our spin networks. We will frequently 
reference the book [|| by Kauffman and Lins; while it explicitly discusses the q- 
deformed version, the formulas we use apply to the classical evaluation as well. 



To set the context, we begin by explaining some elementary algorithms for com- 
puting the lOj symbol. These algorithms all share the feature that they evaluate 
a spin network by choosing bases for the representations labelling the edges, com- 
puting the components of the tensors representing the intertwiners, and computing 
the contraction of the tensors in some way. They make no use of special features 
of these tensors, except for the vanishing property mentioned below. 

The first three methods each have two versions, one which works directly with 
the Spin(4) network (|^), and one that converts it into a five- fold sum over SU(2) 



3. Elementary algorithms 
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Here, U is the spin labeling the new edge introduced by the expansion of the inter- 
twiner at vertex i, and it ranges in integer steps from Li := ma,x{\ji^i—j2,i\, — 
^2,1-2!) to Hi :— min(ji i + J2,i7 ji.i-i +^2, 4-2)1 where the vertex numbers are all to 
be interpreted modulo 5. (If Hi — Li is not a non-negative integer, then the sum 
over li is empty and the lOj symbol is zero. In fact, if vertex i is non-zero, then 
\ji,i-j2,i\, - J2,i-2l, ji,i+j2,i and +j2,z-2 raust all differ by integers.) 

There are at most 0{j) terms in each sum, where j is the average of the ten spins. 
Since the two decagonal networks are the same, one only needs to evaluate one of 
them and square the answer. 

The first elementary method is one we call direct contraction. One simply 
labels each edge in the spin network with a basis vector from the representation 
labelling the edge, and multiplies together the corresponding components of the 
tensors. Then this is summed up over all labellings. In fact, one can restrict to 
a smaller set of labellings: the bases can be chosen so that for each choice of two 
basis vectors on two of the three edges meeting an SU(2) vertex, there is at most 
one choice of basis vector on the third edge giving a non-zero tensor component. 
The Spin(4) vertices also have the property that the bases can be chosen so that 
when three of the basis vectors adjacent to a vertex are specified, the last one is 
determined. 

The second elementary method is staged contraction. In the Spin(4) version 
of this method, one starts with the tensor at vertex 0, contracts with the tensor 
at vertex 1, and then vertex 4, and then vertex 2, and finally vertex 3, again 
taking care to save space and time by using the vanishing properties of the tensors. 
Similarly, one can iteratively contract the tensors in the decagonal SU(2) network. 
At intermediate stages one is storing tensors with a large number of components. 

The third elementary method is Scut. Here one takes a ray from the center 
of (Q) and cuts the three edges it crosses. Then one takes the trace of the operator 
this defines on the three- fold tensor product. In more detail, one sums over basis 
vectors for the factors in this tensor product, computing the effect of the network 
on these basis vectors, and using the vanishing properties. The memory required 
for this method (and the next) is dominated by the memory needed to store the 
tensors themselves. 

The fourth and final elementary method is 2cut. This one only makes sense for 
the decagonal network, since one proceeds by taking a ray from the center of the 
decagon which crosses just two edges, cutting those two edges, and taking the trace 
of the resulting operator, using the vanishing properties. 
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Here is a table which gives an upper bound on the number of operations (addi- 
tions and multiphcations) that these algorithms use, and the amount of memory 
they require, as a function of a typical spin j. The space requirements include the 
space to store the Barrett-Crane tensors. 
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In the last two rows, we represent the entries in a uniform way for cither version 
by writing v = 1 for the Spin(4) version of each algorithm and v = for the SU(2) 
version. Then we let p = j^^^ and m = j5~5v^ This shows how they are related. 
For example, the SU(2) methods always get a factor of in time from the five loops 
coming from expanding the Barrett-Crane vertices. Also, the Spin(4) methods get 
powers of (because dim j <E) j = {2j + 1)^), while the SU(2) methods get powers 
of j (because dimj = (2j -1-1)). 

The 2cut method has the best worst-case behaviour, in both space and time. 

In the next section we present an algorithm which has running time 0{j^) and 
requires 0{p) space, and give variants with running time 0{i^) and 0{i'') and 
which use a constant amount of space. 

The difference between and the running time, j^, for the best of the elementary 
methods is significant. For example, with all spins equal to 20, our algorithm 
runs in under six minutes on a 300 MHz microprocessor. A back of the envelope 
calculation suggests that this would take about 30 years with a algorithm. 

4. Our new algorithms 

In this section we describe a new method for computing lOj symbols. The 
key feature of this method is that it does not proceed by computing the tensor 
components for the intertwiners. Instead it uses recoupling to simplify the network 
to one that can be evaluated directly. Thus this method makes use of special 
properties of the tensors that occur in the lOj symbols. 

There are three versions of this method. We explain one of these in some detail, 
and briefly describe the variants at the appropriate points. We use the notation 
of equation (P, and consider the expansion (^) as a sum of squares of decagonal 
networks. 

The decagonal networks can be deformed into the "ladders" shown below, where 
the vertices at the bottom of each ladder are to be identified with those at the top 
of the same ladder. Any sign introduced by this deformation is produced twice and 
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To simplify these networks further, we recouple the sub-networks consisting of a 
horizontal edge and half of each vertical edge incident to its endpoints, rewriting 
them as sums of sub-networks with the same external edges, using the following 
recoupling formula for SU(2) spin networks |9[ Ch. 7]: 

.6 



E 



b 
d 




The horizontal and vertical networks appearing above are two different ways of 
writing intertwiners from one two-fold tensor product of irreducible representations 
to another; both kinds of networks form bases for the space of intertwiners, and 
the 6j symbols appearing in the formula are defined to be the change-of-basis 
coefficients. 

At first glance, it might look as if this recoupling would introduce sums over ten 
new spins, labeling the ten new vertical edges. However, the total networks that 
result from the recoupling consist of chains of sub-networks shaped like 



(3) 



By Schur's Lemma, such sub- networks will only be non-zero when the incoming 
and outgoing edges have identical spins. So the recoupled networks can be written 
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as a sum over just two new spins, mi and m2. 

fj2,0 ^0 
1^2,4 h 

T. (ft Au.) E {5" 

Zo,--- ,li k=0 mi,m2 

\k J2,3 

I'j2,4 ^4 
1^2,3 ^0 

This sum is over all values such that every vertex in the diagram satisfies the 
Clebsch-Gordan condition, so both mi and m2 will independently range in integer 
steps from maxi(|/i — j2.i-i |) to imiii{li+j2.i-i). Here we use the fact that if the five 
Barrett-Crane vertices are non-zero, then as i varies, the ten quantities \li — j2,i-i| 
and k + j2,i-i aU differ by integers. Indeed, - j2,i-i| = k - J2,i-i = k+ j2,i-i 
modulo integers, and by the paragraph after equation (|^), 

h + hd-l = - J2,i| + i2,i-l = + 324-1 + 32,i = + 32,i+l + 32,i = k+1 + 

modulo integers, where in the third step we use that vertex i -I- 1 is non-zero. 

At this point there is a choice which determines which version of the algorithm 
one obtains. If the sum over mi and m2 is left inside the sum over the Z^, then it can 
be written as the square of a sum over a single m. As described below, each of the 
terms in this sum can be computed with 0{j) operations and a constant amount 
of memory, where j is the average of the ten spins. Thus, this method produces an 
algorithm that runs in ©(j time and takes a constant amount of space. 

In general it turns out to be more efhcient to make the sum over the m's out- 
ermost, in order to reinterpret the sum over the hi as the trace of a matrix prod- 
uct. The range for mi and m2 must encompass all potentially admissible values, 
and the range for each can then be adjusted for the current values of toi and 
m2. The original range for was from Li to Hi (see the paragraph after equa- 
tion (^), so the m's can never be greater than vami^Hi -\- 32,i-i) without violat- 
ing the triangle inequality at one of the vertices. The lower bound is given by 
maxi(mini; (|Zi — j2.i-i| '■ Li < li < Hi)), where the minimum breaks down into 
three cases: if j2,i-i > Hi, it is j2,i-i - Hi] if j2,i-i < Li, it is Li - j2,i-i; and if 
Li < J2,i-i < Hi, it is either or i, depending on whether 2j2,i-i = 2Li (mod 2), 
or not. 

Each of the li is then restricted to take account of the current m values, ranging 
from max(Li, |mi - j2,i-i|, |m2 - j2,i-i\) to min(i?j,mi -|- j2,i-i,m2 +j2A-i)- 

By Schur's Lemma, each sub-network of the form (^ is equal to a multiple of 
the identity, so each of the recoupled ladders, with top and bottom edges joined. 
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is simply a multiple of a loop. KaufFman and Lins |^ give the following formula, 
which can be checked by taking the trace of both sides: 



(4) 



The numerator of the fraction is written 

\2a 




and the denominator is 



(— 1) (2a + 1), the superdimension of the representation. KaufFman and Lins also 
give a formula for the 6j symbol in terms of tetrahedral and 9 networks: 

. -1 Tet[" ; 
a b i[ _ [c d J 

c d 



(5) 

Equations 



.]} 0{a,d,i)e{b,c,i) 
and (|^) allow us to write the following expression for the lOj symbol: 
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(6) 



where 



(7) (AC'"');^' 

The twists implicit in the identification of the top and bottom parts of each network 
introduce signs of (-l)'o+J2,4-r«i ^nd (-i)'o+J2,4-m2 . gj^ce 21q = 2Lo (mod 2), the 
product is (— 1)2(^-0+^2, 4)-mi-ni2 j'j^Q loop values for the spin-mi and spin-rn2 
representations have signs of (—1)2™! and (— 1)^'"^, but since 2mi = 2m2 (mod 2), 
the product of these two signs is always unity. We have made use of the symmetries 
of the tetrahedral networks to put the coefficients in a uniform order for all terms. 
The sum over the k in Equation (||) is the trace of the product of the five matrices 
For each pair of values of mi and m2, these matrices can be computed 
using closed formulas for the tetrahedral and 9 networks given by Kauffman and 
Lins in The formula for the tetrahedral networks involves a sum with 0{j) 
terms, so computing each matrix requires C(j'^) operations]]. The trace of the 
matrix product can also be found in 0{j^) steps. There are two factors of j coming 
from the sums over mi and m2, yielding an overall count of 0{j^) operations. This 
method requires 0{j^) space to store the matrices. 

For some 10-tuples of spins, if all the spins are multiplied by A, the time required 
will scale at a lower power than . Multiplying all the spins by a factor will increase 
the upper and lower bounds of all the sums linearly, but in cases where the two 
bounds are equal, the sum will consist of a single term, regardless of the scaling 
factor. When many of the upper and lower bounds coincide, the first variant of 
the algorithm, with worst case running time ©(j^), in fact becomes faster than the 
0{j^) version. Thus one may wish to use the first variant for certain lOj symbols. 

^Each of the 0{j) terms in the formula for the tetrahedral network contains factorials, which 
themselves require 0{j) operations. However, with some care, the formula can be evaluated with 
a total of 0{j) operations. In practice, wc prccalculate the factorials, using 0{j) space. 
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For large spins, the memory usage can be a problem. For example, with spins of 
around 180, storing each matrix M™^'"^^ requires about 1 gigabyte. In this case, 
one can recalculate the matrix entries as needed, resulting in 0{j^) time and 0{j'^) 
space {0{j^) if factorials are cached). 

The formulas in for the network evaluations are unnormalized. To normal- 
ize all the SU(2) intertwiners according to the convention that any 6 network has 
a value of 1 — which is the convention used in the formula for the Barrett-Crane 
intertwiner — it is simpler to divide the matrix elements by the appropriate 9 net- 
works than to take the existing 9 networks in Equation to be unity and normalize 
the tetrahedral networks. Including this normalization, the matrices become: 

7i^nil,m2\'fc+l 
^'(i2,fc-l, jl,fe) 6'(i2,fe+l, ^/c+l..n,fe+l) 

A subroutine written in C-I-+ that implements this algorithm is available on the 
web along with some sample computations. 

We have not dealt explicitly with the g-deformed case, where the representations 
of SU(2) are replaced with representations of SU(2)q, but it is straightforward to 
adapt each stage of the development above, using the formulas in |^ for the q- 
deformed twist, loop, 9 and tetrahedral networks. 
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